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ABSTRACT 

Aims. Investigation of relationships between dust and gas, and study of the star formation law in M 31. 
Methods. We derive distributions of dust temperature and dust opacity across M 31 at 45" resolution 
using the Spitzer data. With the opacity map and a standard dust model we de-redden the Ha emission 
yielding the first de-reddened Ha map of M 31. We compare the emissions from dust. Ha, HI and H2 by 
means of radial distributions, pixel-to-pixel correlations and wavelet cross-correlations. We calculate 
the star formation rate and star formation efficiency from the de-reddened Ha emission. 
Results. The dust temperature steeply decreases from 30 K near the center to 15 K at large radii. 
The mean dust optical depth at the Ha wavelength along the line of sight is about 0.7. The radial 
decrease of the dust-to-gas ratio is similar to that of the oxygen abundance. Extinction is about linearly 
correlated with the total gas surface density within limited radial intervals. On scales <2kpc, cold 
dust emission is best correlated with that of neutral gas and warm dust emission with that of ionized 
gas. Ha emission is slightly better correlated with emission at 70 ^m than at 24 /im. The star formation 
rate in M 31 is low. In the area 6 kpc < i? < 17 kpc, the total SFR is ~ 0.3 M0yr~^. A linear relationship 
exists between surface densities of SFR and H2. The Kennicutt-Schmidt law between SFR and total 
gas has a power-law index of 1.30±0.05 in the radial range of _R = 7-11 kpc increasing by about 0.3 for 
7? =11-13 kpc. 

Conclusions. The better 70 /xm-Ha than 24/im-Ha correlation plus an excess in the 24 fim/70 fim 
intensity ratio indicates that other sources than dust grains, e.g. of stellar origin, contribute to the 
24 fim emission. The lack of H2 in the central region could be related to the lack of HI and the low 
opacity/high temperature of the dust. Since neither SFR nor SFE is well correlated with the surface 
density of H2 or total gas, other factors than gas density must play an important role in the formation 
of massive stars in M 31. The molecular depletion time scale of 1.1 Gyr indicates that M 31 is about 
three times less efficient in forming young massive stars than M 33. 

Key words, galaxies: individual: M 31 - galaxies: ISM - ISM: dust, extinction - ISM: general -stars: 
formation 



1. Introduction 

Dust, neutral gas and ionized gas are the ma- 
jor components of the interstellar medium (ISM) 
in galaxies. Observations of their properties and 
inter-relationships can give important clues to the 
physics governing star formation. Relationships be- 
tween components in the ISM are to be expected. 
Observations have shown that in the Galaxy dust 
and neutral gas are well mixed. In dense clouds 
of molecular gas mixed with cold dust most of 
the stars are formed. They subsequently heat the 
dust and gas in their surroundings and ionize 
the atomic gas. As the major coolants of the 
ISM are continuum emission and line emission at 
various frequencies, a close comparison of these 
emissions could shed light on spatial and physi- 
cal connections between the emitting components. 
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Present-day IR and radio telescopes have produced 
sensitive high resolution maps of several nearby 
galaxies, which are ideal laboratories to study 
the inter play between the ISM and star forma- 



tion fe.g. iKennicutt et"alll2007t iBigiel et al.ll2008l: 



IVerlev et al.ll2009h 

The spiral galaxy nearest to us, the Andromeda 
Nebula (NGC224), is a highly inchned Sb galaxy of 
low surface brightness. Table 1 lists the positional 
data on M 31. Its proximity and large extent on the 
sky (> 5° X 1°) enable detailed studies of the ISM 
over a large radial range. 

Surveys of M 31 at high angular resolution 
(< 1') are available at many wa velengths. In the 
HI lin e the galaxy was mapped by iBrinks fc Shand 
([1981) at 24" X 36" re solution, the northeast- 
ern half by iBraunI ()1990f ) at 10" resolution and, 
most rec ently, the entire g alaxy with high sensi- 
tivity by iBraun et al.l (|2009D at a resolution of 15". 
iNieten et all (|2006ll made a survey in the -'^^00(1- 
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0) line at a resolution of 23". lDevereux et aD (|1994[ ) 
observed M 31 in the Ha line to obtain the distri- 
bution of the ionized gas. The dust emission from 
M 31 was recently observed by the multiband imag- 
ing photometer Spitzer (MIPS, iRieke et al.ii2004il 
with high sensitivity at 24 fim, 70 fim, and 160 /Ltm 
at resolutions < 40". 

The relationships between gas and dust as well 
as between gas and star formation in M 31 have 
been studie d in the past at resolut i ons o f several 
arcminutes. IWalterbos &: Schwerind (|1987D derived 
a nearly constant dust temperature across M 31 us- 
ing the IRAS 100 fim and 60 /im maps. They also 
found a strong increase in the atomic gas- to-dust 
surface density ratio with increasing radius. This 
increa se was confirmed bv IWalterbos fc Kennicut3 
(|1988D who used optical extinction as dust tracer, 
and by [Nieten et al. (20061 using the ISO map at 
175 /xm ( Haas et al.l 119981 ) Interestinglv. the latter 
authors did not find a radial increase in the molec- 
ular gas-to-dust ratio. 

The dependence of star formation on HI sur- 
face density in M 31 has been studied by a num- 
ber of authors (lEmersonl 1974: Berkhuiisen 1977 : 
Tenies fc HaTIdl 119911: lUnwinI 119801: iNakai fc Sofuel 
1982, 1984) using the number density of HII re- 
gions or of OB stars as star formation tracers. 
They obtained power-law indices between 0.5 and 
2, possibly depending on the region in M 31, the 
star form ation tracer and the angular resolution. 
iBraun et al. (2009) plotted the star formation den- 
sity derived from the brightnesses at 8/im, 24/im 
and UV against the surface densities of molecular 
gas, HI and total gas, but did not fit power laws to 
their data. 

The high-resolution data available for M 31 
show the morphologies of the emission from dust 
and gas components in d etail. We apply a 2-D 
wavelet analysis technique ( Frick et al]|2001f) to the 
MIPS IR data (|Gordon et al.l I2006D and the gas 
(HI, H2, and Ha) maps to study the scale distri- 
bution of emission power and to separate the dif- 
fuse emission components from compact sources. 
We then compare the wavelet-decomposed maps 
at various spatial scales. We also use pixel-to-pixel 
(Pearson) correlations to derive quantitative rela- 
tions not only between different ISM components 
but also between them and the present-day star for- 
mation rate. 



Table 1. Positional data adopted for M 31. 



Following IWalterbos fc Schwerind (|l987t ) and 
iHaas et all (|1998D . we derive the dust temperature 
assuming a emissivity law for the MIPS bands 
at which the emission from the big grains and hence 
the LTE condition is relevant, and present a map of 
the dust color temperature. We also obtain the dis- 
tribution of the optical depth and analyze the gas- 
to-dust surface-density ratio at a resolution of 45" 
(170 pc X 660 pc along the major and minor axis, 
respectively, in the galaxy plane), 9 time s higher 
than before (jWalterbos fc Schwerindll98"7l) . We use 



Position of nucleus 


RA = 00'M2™ 


45.97" 


(J2000) 


DEC = 41° 16 


11.64" 


Position angle of major axis 


37° 




Inclination^ 


75° (0°=face 


on) 


Distance^ 


780±40 kpc^ 





^ Berkhuiisen (1977) and Braun (1991) 
^ Stanek fc Garnavich (1998) 
l'=227±12pc along major axis 



the optical depth map to de- r edden the Ha emission 
observed bv lDevereux et all (|1994D yielding the dis- 
tribution of the absorption-free emission from the 
ionized gas, and use this as an indicator of massive 
star formation. We compare it with the distribu- 
tions of neutral gas to obtain the dependence of 
the star formation rate on gas surface density. 

The paper is organized as follows: The relevant 
data sets are described in Sect. 2. In Sect. 3 we de- 
rive maps of the dust color temperature and optical 
depth, and correct the Ha emission for absorption 
by dust. Radial profiles of the dust and gas emission 
and of the various gas-to-dust ratios are obtained in 
Sect. 4. Sect. 5 is devoted to wavelet decompositions 
and wavelet spectra of the dust and gas distribu- 
tions, and their cross correlations. Complementary, 
we discuss in Sect. 6 classical correlations between 
gas and dust. In Sect. 7 the dependence of the star 
formation rate on the gas surface density is pre- 
sented. Finally, in Sect. 8 we summarize our results. 



2. Data 

Table 2 summarizes the data used in this work. 
M 31 was mapped in IR (at 24 /im, 70 /im, and 
160 fim) by MIPS i n August 2004 cover ing a region 
of about 1° X 3° (jCordon et al.ll2006D . The basic 
data reduction and mosaicing was performed using 
the MIPS instru ment team Dat a Analysis Tool ver- 
sions 2.90 (Gor don et al.ll2005D . 

M 31 was observe d in the ^^00 ( 11-0) line with 
the IRAM telescope bv lNieten eraD (12006D at a res- 
olution of 23". They derived the distribution of the 
molecular gas using a constant conversion factor of 
Xco = 1-9 X 10^" mol. K^^ km~^ s. The galaxy was 
observed in the 2 1-cm HI line with the W esterbork 
interferometer bv lBrinks fc Shan^ ()1984l ) at a res- 
olution of 24" X 36". The HI survey has been cor- 
rected for missin g spacings. The Ha observations of 
iDevereux et al.l (|1994) were carried out on the Case 
Western Burell-Schmidt telescope at the Kitt Peak 
National Observatory, providing a 2° x 2° field of 
view. 

Although the resolution of 40" of the 160 /im 
image is the lowest of the data listed in Table 2, 
we smoothed all maps to a Gaussian beam with a 
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half-power width of 45" for a comparison with ra dio 
continuum data at 20 cm (jHoernes et al.lll998[ ) m 
a forthcoming study (Tabatabaei et al. in prep.). As 
the point spread function (PSF) of the MIPS data 
is not Gaussian, we convolved the MIPS images us- 
ing custom kernels created with Fast Fourier trans- 
forms to account for the detailed structure of the 
P SFs. Details of the k ernel creation can be found 
in iGordon et all (l20Q7h . 

After convolution, the maps were transformed 
to the same grid of 15" width with the reference 
coordinates and position angle of the major axis 
given in Table 1. Finally, they were cut to a com- 
mon extent of 110' x 38. '5, for which most data 
sets are complete. The field is not centred on the 
nucleus of M 31, but extends to 56. '25 along the 
northern major axis (corresponding to a radius of 
R = 12.8 kpc) and to 53. '75 along the southern ma- 
jor axis {R = 12.2 kpc). The II2 map of Nieten et al. 
(2006) extends to 48. '5 along the southern major 
axis {R— 11.0 kpc). With an extent of 19. '25 along 
the minor axis in both directions, the field covers 
radii of i? < 16.9 kpc in the plane of M 31. Hence, 
radial profiles derived by averageing the data in cir- 
cular rings in the plane of the galaxy (equivalent to 
elliptical rings in the plane of the sky) are incom- 
plete at i? > 12 kpc because of missing data near 
the major axis. 

3. Dust temperature and opacity 

IWalterbos fc Schwerind (jl987D extensively studied 
the distributions of the dust temperature and opac- 
ity in M 31 using the IRAS data at 60 fim and 
100 /im smoothed to a resolution of 4. '3 x 6. '9. 
Assuming a emissivity law, they found a 

remarkably constant dust temperature (21-22 K) 
across the disk between 2 kpc and 15 kpc radius. 
Using this temperature, they obtained the opacity 
distribution at 100 fim. Below we apply a similar 
method to the 70 /im and 160 /xm MIPS maps to 
derive the distributions of dust temperature and 
optical depth at the Ha wavelength at higher reso- 
lution and sensitivity. 

3.1. Dust temperature 

We derived the color temperature of the dust, T^, 
between 70 /im and 160 ^m assuming a emis- 
sivity law that should be appropriate fo r interstel- 
lar grains emitting at the se wavelengths (jAndriessd 
I1974I : iDraine fc Leel[T98l . The resulting dust tem- 
perature map (Fig. 1) and a histogram of the tem- 
peratures (Fig. [2^) show that Td varies between 
15.6 ± 0.8 K and 30.8 ± 0.3 K. The mean value 
of 18.7 ± 1.4 K (sta ndard deviation) is low e r than 
that obtained by IWalterbos fc SchwerinS ()1987t ) 
and close t o the ISO measurements (16±2K) of 
iHaas et all (|1998[ ). Figure 1 shows that dust of 
~18K exists all over M 31. Warmer dust with 



Td > 20 K dominates in star forming regions and in 
an extended area around the center of the galaxy, 
while cooler dust dominates in interarm regions. 

Figure |2Jd shows the dust temperature averaged 
in rings of 0.2 kpc width in the plane of M 31 against 
radius RQ. On both sides of the center the dust 
temperature falls very fast from about 30 K near 
the nucleus to 19 K at i? ~ 4.5 kpc. To the outer 
parts of the galaxy, it then stays within a small 
range of about 17K-19K in the north and 16 K- 
19 K in the south. This indicates different radiation 
characteristics between the inner 4 kpc and beyond. 
In the ring of bright emission, the so called 
'10 kpc ring', the temperature is clearly enhanced, 
especially in t he northern half. Thus, in cont rast to 
the finding of lWalterbos fc Schwerinel ([1983), is 
not constant in the range i? = 2-15 kpc but varies 
between 22.5 ± 0.5 K and 17.2 ± 0.7 K. 

It is interesting that the mean dust tempera- 
ture obtained between 70 fim an d 160 nm. is about 
3K low er in M 31 than in M 33 (jTabatabaei et all 
l2007b[ ). The emission from cold dust in M 31 is 
stronger than in M 33, which can also be in- 
ferred from the total emission spectra based on 
IRAS and ISO observations (see iHaas et al]ll998l : 
iHippelein et~aL| [2003). 

3.2. Dust opacity distribution 

The total dust optical depth along the line of 
sight T160 was obtained from the dust intensity 
at 160 jCtm and th e derive d temperature. Following 
iTabatabaei et al.l (|2007bD . T160 was converted into 
the dust optical depth at the wavelength of the 
Ha line, tho , by multiplying it by the ratio of the 
dust extinction coefficient per unit mass at the cor- 
responding w avelengths. T ug — 2200 ngo (see e.g. 
Figure 12.8 of lKriige]|l2003D . Figure 3 shows the dis- 
tribution of THa across the disk of M 31 at an an- 
gular resolution of 45". Regions with considerable 
dust opacity {tuo > 0.6) follow the spiral arms, 
even the inner arms which are either weak or not 
detected in Ha emission. The high opacity clumps 
(thq > 1-5), however, only occur in the arms at 
~5kpc and in the 'lOkpc ring'. On average, the 
optical depth is largest in the TO kpc ring'. Hence, 
in this ring, dust has the highest density like the 
atomic gas ()Brinks fc Shanel ll984D . Fig. Hi shows 
the histogram of THa, indicating a most probable 
value of 0.5 and a mean value of 0.7 ±0.4 across 
the galaxy. This agrees with the value of thq = 
0.5 ± 0.4 that follows from the mean total extinc- 
tion obtained by Barmby et al (2000) towards 314 
globular clusters. 



Because the spiral structure is different in the north- 
ern and southern half (northeast and southwest of the 
minor axis, i.e. left and right of the minor axis, respec- 
tively), we present all radial profiles for each half sepa- 
rately. 
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Table 2. M 31 data used in this study. 



Wavelength 


Resolution 


Telescope 


Ref. 


160 fim 


40" 


Spitzer 


Gordon et al. (2006) 


70 /im 


18" 


Spitzer 


Gordon et al. (20^) 


24 /im 


6" 


Spitzer 


Gordon et al. (2006) 


2.6nini^2CO(l-0) 


23" 


IRAM 30-ni 


Nicten et al, (2006) 


21 cm HI 


24" X 30" 


WSRT 


Brinks & Shane (1984) 


6570AHa 


2" (pixel size) 


KPNO 


Dcvereux et al. ('1994) 




11 50 
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Fig. 1. Dust temperature in M 31 obtained from the ratio lyo/^m/Iieo^im based on the Spitzer MIPS data. 
Only pixels with intensity above 3 x noise level were used. The angular resolution of 45" is shown in the 
lower right-hand corner of the map. The cross indicates the location of the center. The bar at the top 
gives the dust temperature in Kelvin. 



The variation of the mean dust optical depth 
TRa with galactocentric radius is shown in Fig. HJd. 
In the north, tho peaks not only in the '10 kpc ring' 
(with two maxima at R =9.9 and 10.9 kpc) but 
also near 5 kpc (with two maxima at R =4.3 and 
5.9 kpc). Beyond 11 kpc thq drops with an expo- 
nential scale length of 2.48±0.07kpc in the north 
and 5.06±0.22kpc in the south (Sect. 4, eq. 1). 

In Fig. [S] we compare the radial variation of 
tho! for the total area with earlier determina- 
tions. The various estimates ag r ee we ll given the 
large uncertainties. IXu &: Heloul ()1996D derived ry 
from high-resolution IRAS data using a dust heat- 
ing/cooling model and a sandwich configuration of 
dust and stars. Although they left out the dis- 
crete sources, their values may be too high be- 
cause they did not include inter-arm regions in 
their study. iMontalto et al.l (pOOg) calculated the 
extinction Apuv from the total infrared TIR-to- 
FUV intensity ratio and a sandwich model for stars 
and dust. They note that at i? < 8 kpc the geome- 
try of M 31 may differ from the sandwich model 
due to the stars in the bulge, making the inner 
points less reliable. This would also affect the re- 



sults of lXu fc Heloul (|1996[1 at i? < 8 kpc. The TIR- 
to-FUV ratio is applicable if the dust is mainly 
heated by young stars, but in M 31 a bout 70% of 
the c old dust is heated by the ISRF (|Xu fc Heloul 
|1996[) . Therefore, A puy is overestima ted, as was 
al so found for M 33 (Vcrl ev et al1 l200^. The curve 
of lTemoel et"all (JOIQ) closely agrees with our data. 
They derived ry from MIPS data and a star-dust 
model. Their smooth curve underestimates ry in 
the brightest regions by about 0.1 and may overes- 
timate Ty in regions of low brightness. 

The opacity map in Figure 3 can be used to 
correct the Ha emission for the extinction by dust. 
In general, extinction depends on the relative dis- 
tribution of emitting regions and dust along the 
line of sight and changes with the geometry (e.g. 
well mixed diffuse medium or shell-like in HII re- 
gions IWitt k Gordon 200Q ). In this study, individ- 
ual HII regions are rarely resolved and the geometry 
is close to a mixed diffuse medium. Furthermore, 
there is no information about the relative position 
of emittors and abs orbers along the lineof sight. 
For the Milky Way, iDickinson eFaP (|2003t) found 
indications of a non-uniform mixing by comparing 
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Dust temperature (K) R (Kpc) 

Fig. 2. a) Histogram of tlie dust temperature sliown in Fig. 1. b) Distribution of the dust temperature 
in rings of 0.2 kpc in the galactic plane in the northern and southern halves of M 31. 



the z-distribution of atomic gas and dust. They 
adopted one third of the total dust optical depth 
as the effective extinction as a first-order ap proach. 
This is also in agreement with lKruigel (l2009l) taking 
scatte ring into account. Moreover, Magnier et al.l 
(|1997D found that on average the extinction comes 
from dust associated with only one-third of aver- 
age N(HI) in their study of OB associations along 
the eastern spiral arm regions of M 31. Therefore, 
we use an effective optical depth Tcs = 0.33 x thq 
in this paper. The attenuation factor for the Ha 
intensity then is e""^"*' and we derive the intrin- 
sic Ha intensity Iq from the observed Ha intensity 



Integration of the Ha map out to a ra- 



dius of 16 kpc yields a ratio of corrected-to-observed 
total Ha flux density of 1.29, thus about 30% of the 
total Ha emission is obscured by dust within M 31. 
The corrected Ha map is shown in Fig. llOb ,. 

Near the center {R < 1 kpc) , Tcs varies be- 
tween 0.03 and 0.13, corresponding to an extinc- 
tion of A{Ha) = 1.086 x Tcff ~ 0.03-0.14 mag. 
At larger radii, the mean extinction increases, par- 
ticularly in dense clouds and star forming regions, 
reaching a maximum of yl(Ha) ~ 1.2 mag at the 
densest dust cloud in the south-east of the '10 kpc 
ring' (RA^OO'' 41™ 05.10^ and DEC = -^40° 38' 
17.73"). The range of extinction values agrees with 
that derived from the op t ical s tudy of dust lanes by 
IWalterbos fc KennicuttI (Il988f) and the photomet- 
ric studv of lWilliamsl (|2003D . 



4. Radial distributions of dust and gas 
emission 

4.1. Radial profiles 

In this section, we present the mean surface bright- 
ness along the line of sight of dust and gas com- 
ponents as a function of galactocentric radius R. 
The surface brightnesses are averaged in 200 pc- 



wide circular rings about the nucleus in the plane 
of M 31. This is equivalent to averaging in elliptical 
rings of 53" width in the plane of the sky. 

For simplicity we used a constant inclination an- 
gle of 75° at all radii, appropriate for the emission 
at i? > 30' (6.8 kpc), although in H a and HI the in- 
ner r e gions are seen more face-on (ICiarduUo et all 
119881: iBraunI 119911: IChemin et all l2009h . However, 
using i = 75° for R < 30' instead of i = 68° (the 
area- weighted mean of the incli nations for the inter- 
vali? = 1.9-6.8kDC given bv [Chemin et al.ll2009t) 
does not change our results. The smaller inclination 
shifts the radial positions of the inner arms about 
0.5 kpc inwards but the general shape of the pro- 
files remains the same, and as all profiles change in 
a similar way their inter-comparison is not affected. 
Furthermore, the results of the classical correlations 
for R < 30' presented in Sect. 6 are the same within 
the errors for i = 68° and i = 75°. 

Figure [6] shows the mean IR intensities and the 
gas surface densities versus the galactocentric ra- 
dius R for the northern and southern halves of 
M 31. The radial profiles of the IR emission at 
24 /im and 70 fim are similar. The 160 /im emis- 
sion, representing the colder dust emission, how- 
ever, shows a generally flatter radial distribution 
than the 24 /im and 70 /im emission. In particu- 
lar, the fast decrease of the 24 /im and 70 fim pro- 
files from the center to i? ~ 2 kpc doe s not occur 
at 160 /im. This is in agreement with iHaas et "all 
(jl998f ) who concluded from their ISO 175 /im map 
and IRAS data that the dust near the center is rel- 
atively warm. The fast central decrease of warmer 
dust emission may be attributed to a decrease in 
the UV radiation field outside the nucleus, as a 
similar trend is seen in the GAL EX UV profiles 
presented bv iThilker et al.l (|2005f ). At all three IR 
wavelengths the arms are visible, even the weak in- 
ner arms. The bright arms forming the '10 kpc ring'. 
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Fig. 3. Distribution of the dust optical depth at Ha wavelength in M 31. The bar at the top shows th„. 
The angular resolution of 45" is shown in the lower right-hand corner of the map. The cross indicates the 
location of the center. Overlayed are contours of molecular gas column density N(H2) with levels of 250 
and 800 xlO^^ mol. cm^^. Note that maxima in not always coincide with maxima in N(H2). 




Fig. 4. a) Histogram of the dust optical depth shown in Fig. 3, 6) radial distribution of the mean optical 
depth at the Ha wavelength in rings of width of 0.2 kpc in the galactic plane in the north and south of 
M 31. The errors are smaller than the size of the symbols. 



are pronounced in the north and followed by an ex- 
ponential decrease toward larger radii. 

Although the general trend of the warm dust 
surface brightness (at 24 /xm and 70 /xm) resembles 
more that of Ha than those of the neutral gas pro- 
files (Fig. El lower panels), small variations (e.g. in 
the inner 5 kpc and for R— 11- 12 kpc in the south) 
follow those in the total gas distribution due to vari- 
ations in the molecular gas. Beyond about 5 kpc, 
the radial profile of the cold dust (160 /im) is simi- 
lar to that of the molecular gas, but with smoother 
variations. The minimum between 5 kpc and 10 kpc 
radius at 24 /xm and 70 /xm is less deep at 160 /im 
and is missing in the HI profile. 



We obtained radial scale lengths between the 
maximum in the '10 kpc ring' and R=14.9 kpc for 
the northern (Zn) and southern (Is) halves of M 31 
separately as well as for the total area (/). We fit 
an exponential function of the form 

I{R) = Ioexp{~R/l), (1) 

where Iq is the intensity at i? = 10.9 kpc for the 
total area and in the north, and R — 8.9 kpc in 
the south. The resulting scale lengths are listed in 
Table O 

In each half of M 31, the scale lengths of 
the warm dust emission are smaller than that of 
the cold dust. This confirms that the warm dust 
is mainly heated by the UV photons from the 
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Fig. 6. Top: Radial profiles of the Spitzer IR emission from the northern (left) and the southern (right) 
halves of M 31. Bottom: Radial profiles of the surface densities of the atomic, molecular and total neutral 
gas together with that of the ionized gas (de-reddened Ha) for the northern (lejf) and southern {right) 
halves of M 31. The units are 10^^ atoms cm~^ for HI and HI-I-2H2, 10^* molecules cm~^ for H2, and 
10^*^ ergs~^ cm~^ sr~^ for Ha profiles. The profiles show intensities along the line of sight averaged in 
circular rings of 0.2 kpc width in the plane of M 31 against the galactocentric radius. The errors are 
smaller than 5% for all profiles, only for H2 they increase from 10% to 25% at i? > 12.3 kpc and in the 
inner arms for R < 4.5 kpc. 



star forming regions in the '10 kpc ring' and the 
cold dust mainly by th e interstellar radia tion field 
(ISRF) from old stars (iXu fc Heloul[T996h . 

The scale lengths of the 24 /j,m and 70 iim emis- 
sion are nearly the same and the 24 /im-to-70 /xm 
intensity ratio (Fig. [T]) hardly varies between R ~ 
10 kpc and R ~ 15 kpc. This indicates a similar dis- 
tribution of their origins. Assuming that the main 
source of the 24 //m emission is very small dust 
grains and of the 70 /xm and 160 /xm emission is 
big gr ains, as argued by IWalterbos fc Schwerind 
(|1987[ ). a constant intensity ratio of the 24/xm-to- 
70 /xm and 24/im-to-160 /im emission suggests that 
the very small and big grains are well mixed in 
the interstellar medium. Other possible origins of 
the 24 /im emission are stars with dust shells like 
evolve d AGB stars or C arbon stars. Using the IRAS 
data, (jSoifer et al.lllQSSl attributed the 25 ^m emis- 
sion from the bulge (central 8') of M 31 to cir- 



cumstellar dust e mission from late-type stars. In 
the disk of M 31, IWalterbos fc Schwerind (|1987D 
found no direct evidence for a contribution from 
stars with dust shells ( contrary to that in the Milky 
Way, ICox et al.lfl98fih . 

The higher resolution and sensitivity of the 
MIPS IR intensity ratios (Figl?]), however, provide 
more information. Although at > 3 kpc the vari- 
ations in the IR intensity ratios are not large, their 
radial behavior is not the same. For instance, the 
24-to-70 /im intensity ratio peaks between 5 kpc and 
10 kpc radius, whereas the 70-to-160/im intensity 
ratio peaks in the '10 kpc ring'. The latter can be 
explained by the higher temperature of the dust 
heated by OB associations in the '10 kpc ring'. 
The fact that the 24-to-70 /im intensity ratio is not 
enhanced in the '10 kpc ring' (and in the central 
region) shows the invalidity of this ratio for tem- 
perature determination due to the important con- 
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RADIUS [kpc] 

Fig. 5. Radial variation of the (total) mean op- 
tical depth in Ha along the line of sight for the 
full area (north+south). Plusses: our data averaged 
in 0.2 kpc- wide rings in the plane of M 31 using 
i = 75°; stars: same but with i = 77.5° for compar i- 
son with other work; triangles: IXu fc Heloul ()1996[ ). 
averages in 2 kpc- wi de rings with i = 77° (sc aled to 
D = 780kpc); circles: iMontalto et alJ (|2009l ). aver- 
ages in 2 kpc - wide rings with i = 77.6°; solid line: 
iTemoel et all (|2010t) . semi-major axis cut through 
model with i = 77.5°. The errors in our data and in 
the curve of iTempel et alJ ()2010D are about 10% of 
the mean values. 



Table 3. Exponential scale lengths of dust and gas 
emissions from M 31. The scale lengths were cal- 
culated from i?= 10.9 kpc to i?= 14.9 kpc for the 
whole galaxy (l) and the northern half (Zn) and 
from i?=8.9kpc to i?=14.9kpc for the southern 
half (Is)- The scale length of rna is also shown for 
comparison. 





In (kpc) 


Is (kpc) 


I (kpc) 


IR leO^im 


1.86 ±0.06 


3.87 ±0.18 


2.29 ±0.11 


IR 70^m 


1.38 ±0.08 


3.07±0.17 


1.66 ±0.09 


IR 24^m 


1.38 ±0.07 


3.44 ±0.21 


1.57 ±0.06 


H2 


1.14 ±0.07 


2.40 ±0.22 


1.27±0.11 


HI 


4.54 ±0.24 


15.90 ±1.70 


5.90 ±0.20 


HI-F^2H2 


3.80 ±0.12 


9.00 ±0.69 


4.73 ±0.21 


Hq 


1.08 ±0.03 


2.92 ±0.10 


1.39 ±0.04 




2.48±0.07 


5.06±0.22 


3.17±0.15 



tribution from the very small grains. On the other 
hand, the enhancement of the 24-to-70 fim in re- 
gions where there is no strong radiation field (be- 
tween the arms) reveals possibly different origins 
of the 24 /im and 70 /im emission. The stellar ori- 
gin, e.g. photosphere of cool stars or dust shell of 
the evolved stars, may provide the enhancement 
of the 24-to-70 /im intensity ratio in the inter-arm 
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Fig. 7. Ratio of the MIPS IR intensities against 
galactocentric radius in M 31. 



region. In M 33, IVerlev et all (|2009[ ) attributed a 
similar enhancement of the diffuse 24 /zm emission 
to dusty circumstellar shells of unresolved, evolved 
AGB stars. For M 31, this needs to be quantified 
through a more detailed study and modelling of the 
spectral energy distribution, which is beyond the 
scope of this paper. 

4.2. Gas-to-dust ratio 

The gas-to-dust mass ratio and its variation across 
the galaxy can provide i nformation about the 
metallicity distribution (e.g. IViallefond et al.lll982[ ) 
and hence about the evolutionary history of the 
galaxy. The relative amount of dust and gas is ex- 
pected to be cor related with the ab undance of the 
heavy elements (jPraine et al.ll2007ll . 

A number of authors has studied the gas-to-dust 
ratio in M 31 by comparing HI column densities and 
optica l or UV extinct ion (Waltcrbos & Kcnnicutlj 
1988; Baiaia fc Gereelvi 119771: iXu fc Hclou 199e' 



Nedialkov et al.l 120001: ISavcheva fcTassev .20d' 



All authors found an incr ease of the atomic gas- 
to-dus t ratio with radius. 'Walte rbos &: Schweringj 
(|l987t) derived the HI gas-to-dust ratio using dust 
optical depth from IRAS 60 /im and 100 fim data. 
They found a radial gradient th at is 4-5 time s large r 
than the abundance gradient of I Blair et al.l (|l982l) . 
After adding the molecular and atomic gas column 
densities, Nieten et al. (2006) obtained a strong ra- 
dial increase in the total gas-to- 1 75 /im intensity 
resulting from the increase in the atomic gas-to- 
175 iim intensity. As the dust optical depth is a 
better measure for the dust column density than 
the temperature-dependent dust emission, we re- 
investigated the gas-to-dust ratio in M 31 taking ad- 
vantage of the high resolution of the Spitzer MIPS 
data. 

We calculated the radial profiles of the three 
gas-to-dust ratios from the mean column densities 
of N(HI), N(2H2), N(HI+2H2) and ma in circular 
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Fig. 8. Radial profiles of the gas-to-dust ratios in 
M 31, the northern half and the southern half. 
Top: N(HI)/THa, middle: N(2H2)/thq , bottom: 
N(HI+2H2)/rHQ. In the middle panel, the north- 
ern profile is shifted by 300 units for clarity. The 
errors are smaller than 5% everywhere only for 
N(2H2)/thq they increase from 10% to 25% at 
R > 12.3 kpc and in the inner arms for R < 4.5 kpc. 



rings of 0.2 kpc width in the plane of the galaxy. 
Figure [S] (upper panel) shows that the atomic gas- 
to-dust ratio increases exponentially with radius by 
more than a factor of 10 from about 0.6 x 10^"'^ 
atcm~^ at the center to about 6.5 x 10^^ atcm~^ 
at i?=15kpc. The increase is surprisingly smooth 
and, at least up to R= 13 kpc, nearly the same for 
the northern and southern half, indicating little 
variation between arm and inter-arm regions and 
within the arms. In contrast, the molecular gas- 
to-dust ratio (Fig. [51 middle panel) does not in- 



crease systematically with radius but shows clear 
enhancements of a factor 2-3 in the spiral arms 
and the '10 kpc ring'. The minima in the inter-arm 
regions are due to a stronger decrease in N(2H2) 
than in rna- Figure 3 shows that along the arms 
N(2Il2)/rHQ also varies significantly because max- 
ima in II2 emission and thq are often not coinci- 
dent. The variations in N(2H2)/rHa are visible in 
the profile of the total gas-to-dust ratio (Fig. [51 bot- 
tom panel) as weak enhancements at the positions 
of the arms near i? = 6 kpc and R — 8 — 12 kpc. As 
the atomic gas is the dominant gas phase in M 31, 
dust mixed with HI gas largely determines the opti- 
cal depth. Inspection of the distribution of the total 
gas-to-dust ratio across M 31 (not shown) reveals 
small-scale variations along the arms of typically a 
factor of 2. 

We conclude that the radial increase in the total 
gas-to-dust ratio of more than a factor 10 between 
the center and R = 15 kpc is entirely due to that 
of the atomic gas-to-dust ratio, whereas the molec- 
ular gas-to-dust ratio is only incre ased in the arms . 
This confirms the conclusion of Nieten et al.l (|2006f) 
based on the same gas data and the 175 /xm inten- 
sity. 

At which radius in M 31 would the gas-to-dust 
ratio observed in t he sol ar neighborhood occur? 
iBohlin et al.l (|1978[) and iDinlas fc Saw3 (|1994ll 
derived N(HI)/E(B-V)=4.8 x lO^^ atcm-^mag-^ 
and (4.9 ± 0.3) x 10^^ atcm~^mag~\ respec- 
tively, using the extinction towards large sam- 
ples of stars to determine the color excess E(B- 
V). Since E(B-V) = Ay/Rv, where the visual ex- 
tinction Ay = 1.234 THa mag (e.g. iKriigell [2003[l 
and the total/s elective extinction Ry =2.8ib0.3 
in M 31 (Walterbos fc Kennicut"3 Il988^ ■ we have 
E(B-V)=0.44THa mag~^. Hence, a value of 
N(HI)/E(B-V)=4.9 X lO^i atcm-^mag-^ corre- 
sponds to N(HI)/rHa = 2.2 x 10^^ atcm'^, which 
occurs in M 31 near R =8.5 kpc (Fig. [51 top 
panel), just in the bright emission ring. The to- 
tal gas-to-dust ratio near the sun of 5.8 x 10^^ 
atcm~^mag~^ (jBohlin et al.l [T9781 ) corresponding 
to N(gas)/rHa= 2.6 x 10^^ at cm~^ occurs at nearly 
the same radius (Fig. [51 bottom panel). Thus 
the gas-to-dust ratio near the sun is similar to 
that in the '10 kpc ring' in M 31, in agree- 
ment with earlier studies (Ivan Cxendered 119731 : 
IWalterbos fc Schwering 1983). 

In contrast to Fig. [H we present in Fig. 9 the 
radial profiles of the dust-to-gas ratios, here for the 
total area in M 31. The two lower curves closely fol- 
low exponentials with scale lengths of 6.1±0.2kpc 
and 7.4±0.2kpc for THa/N(HI) and THa/N(gas), 
respectively, between R= 5 kpc and R= 15 kpc (see 
Table 4) . For nearly the same radial ra nge (i?=3- 
15 kpc), IWalterbos fc Schweringj (|1987t ) derived a 
scale length of T -ioomn/N(HI)~4kpc from data near 
the major axis. IWalterbos fc KennicuttI (|1988D ob- 
tained a scale length of Ab/N(HI)~ 9kpc for the 
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Table 4. Exponenti al scale leng ths L and r adial gradients of du s t-to-g as ratios and the abundance [0/H] 
combined data from lBlair et alj (|1982f ) and iDennefeld fc KunthI ()198H ). where Ratio(i?) = C. exp(-i?/L). 
N(HI) and N(gas) are in lO^^atcm"^. Errors are standard deviations. 



Ratio 


R 




C 


L 


Gradient 




(kpc) 






(kpc) 


(dex/kpc) 


rH<./N(HI) 


0-16 


1.95 


±0.08 


6.2±0.1 


0.070 ±0.001 




5-15 


1.97 


±0.08 


6.1±0.2 


0.071 ±0.001 


-rHc«/N(gas) 


0-16 


1.35 


±0.05 


7.3 ±0.2 


0.059 ±0.002 




5-15 


1.27 


±0.06 


7.4 ±0.2 


0.059 ±0.002 


[O/H] X IQ-'' 


5-15 


2.95 


±0.35 


9.7±2.6 


0.045 ±0.012 



in ner and outer dust lanes, and the tv/N(HI) ratio 
of IXu fc Heloul (|l996) for diffuse spiral arm regions 
also indicates a scale length of about 4 kpc (all scale 
lengths were scaled to D= 780 kpc) . Since our scale 
lengths are not restricted to specific areas, our re- 
sults are more representative for the mean dust-to- 
gas ratios in the disk of M 31. 

As dust consists of heavy elements and both 
dust and heavy elements are found in star forma- 
tion regions, the radial variations in the dust-to- 
gas ratio and the metal abundan ce are expected 
to be similar (e.g. (jHirashita 1999i; iHirashita et al.l 
I2OO2I) . This has indee d been observ ed in sev- 
eral nearby galaxies (|Issa et al.l |1990() . In M 31 
the variation in the metallicity with radius is 
not well established. Measurements of the ele- 
ment abundance strongly depend on the empiri- 
cal rn e thod and ca. li bratio n applied (j Trundle et al.l 
L2OO2I). i Pagel et al.l (Il979h showed that the ([Oil] 
±[OIII])/H/? ratio (the so called 'R23') is a good 
probe of oxygen abundance and radial trends 
of this ratio ha ve b een studied in many nearby 
galaxies (e.g. | Pag cl & Edmun dj 1981 ; lEvarisI 
19861: iHe nrv fc Howard,1995; .Garnett et al.l^l997^ ■ 
Blair et al. ( 1982f) and IDennefeld fc KunthI (|l98lh 
derived R23 for HII regions in M 31. We com- 
bined their results and derived a scale length of 
log[0/H] of 9.7±2.6kpc corresponding to a gradi- 
ent of 0.045 ± 0.0 12 dex/kpc. Compari ng four differ- 
ent calibrations, iTrundle et al.l (|2002D derived gra- 
dients o f 0.027- 0.0 13 dex/kpc using the 11 HII re- 
gions o fH air et al.| (fl982 ). with Pagel's calibration 
giving 0.017 ±0.001 dex/kpc. Since our value of the 
[0/H] gradient is based on 19 HII regi ons, we expect 
it to be more reliable than that of ITrundle et al.l 
(iOOl). 

Table 4 shows that the radial gradient in 
THQ/N(gas) best matches the metallicity gradient. 
In view of the large uncertainties, the gradients in 
the dust-to-gas surface-density ratio and the oxygen 
abundance in M 31 may indeed be comparable. A 
much larger sample of abundance measurements of 
HII regions is needed to verify this similarity. Our 
result agrees with the approximately linear trend 
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Fig. 9. Dust-to-gas ratios as function of galac- 
tocentric radius for M 31, calculated from 
the radial profiles of TRa, N(HI), N(2H2) and 
N(HI±2H2)= N(gas). 



between gradients in dust-t o-gas ratios and [O /H] 
in nearby galaxies noted bv llssa et al.l (|l990f ). 



5. Wavelet analysis of dust and gas emission 

To investigate the physical properties of different 
phases of the interstellar medium as a function of 
the size of emitting regions, wavelet transforma- 
tion is an ideal tool, yje use the Pet Hat wav elet 
(see iFrick et al.ll200lt iTabatabaei et al.ll2007a^ to 
decompose the emissions of IR, HI, H2, HI±2H2, 
and de-reddened Ha into 10 spatial scales starting 
at 0.4 kpc (about twice the resolution). The cen- 
tral 2 kpc was subtracted from all images before 
the wavelet transformation to prevent a strong in- 
fluence of the nucleus on the results. As an example, 
we show the extinction-corrected Ha map and the 
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Ha emission for 3 different scales in Fig. \TU[ On 
the scale of 0.4 kpc, the distribution of HII com- 
plexes and large HII regions is borne out. The scale 
of 1.6 kpc (the typical width of spiral arms) shows 
connected HII complexes along the arms, and on 
the scale of 4 kpc we see the extended emission from 
the '10 kpc ring'. 

5.1. Wavelet spectra 

The wavelet spectrum, M(a), represents the distri- 
bution of the emitting power as function of the scale 
a. The wavelet spectrum will smoothly increase to- 
wards larger scales if most of the emission is coming 
from diffuse structures forming the largest scales, 
here up to 25 kpc. On the other hand, the spec- 
trum will decrease with increasing scale if compact 
structures are the dominant source of emission. The 
spectra of the IR and gas emission are shown in 
Fig.ini 

All IR and gas spectra are intermediate between 
the two cases described above. Only the spectra of 
the HI gas and the 160 ^m emission generally in- 
crease with scale indicating the importance of dif- 
fuse HI and cold dust emission. In addition, the HI 
spectrum exhibits a dominant scale at a ~ 4 kpc 
corresponding to the width of the '10 kpc ring', 
where strong diffuse emission occurs in interarm re- 
gions. The large width of the HI 'ring' is also visible 
in the radial profiles in Fig. [51 The dominant scale 
of the emission from warm dust, molecular gas and 
Ha is near 1 kpc, where complexes of giant molec- 
ular clouds and star forming regions show up. The 
IR spectra at 24 /j,m and 70 /im on scales a < 6 kpc 
look most similar indicating that the star forming 
regions are the main heating sources at both wave- 
lengths. On the other hand, the effect of the ISRF 
heating the cold dust is well indicated in the 160 fim 
spectrum where a general increase towards larger 
scales is found. All spectra, apart from that of HI, 
show a minimum near a = 6 kpc corresponding to 
the large, weak interarm region inside the '10 kpc 
ring'. The spectrum of Ha is most similar to that of 
70 /im, which may explain why the Ha emission cor- 
relates better with 70 /im emission than with that 
at 24 /im (see Sect. 6.2 and Table 6). 

The spectrum of the Ha emission is flat on small 
scales up to 1.6 kpc, the width of the spiral arms in 
the Ha map. This is understandable as the emission 
from very compact HII regions is unresolved at our 
resolution and not many large HII complexes exist 
especially in the south (see the decomposed map in 
Fig.Iinb of a==0.4kpc). 

5.2. Wavelet cross-correlations 

We derive the cross-correlation c oefficients, rui(a), 
for dif ferent scales following iTabatabaei et all 
(|2007a[ ). The correlation coefficients are plotted in 
terms of scale in Fig. [T^l They show that IR emis- 



sion correlates with the emission from different gas 
phases on most scales. In all cases, emission from 
structures on scales larger than 10 kpc are best cor- 
related. This corresponds to scales of the diameter 
of the '10 kpc ring' and the over-all structure of the 
galaxy. On medium scales, the weakest correlation 
occurs between HI and dust emission on a = 6 kpc. 
This scale includes areas of significant diffuse HI 
emission where the dust emission is weak interior 
to the 'lOkpc ring' (compare also Fig. [TT|) . On the 
smallest scale of 0.4 kpc, the cold dust emission is 
best correlated with that of the total neutral gas, 
while the warm dust emission at 70 /im is best cor- 
related with the ionized gas emission. Note that on 
this scale, the 24 /im and 70 /im (warm dust) emis- 
sions hardly correlate with HI (^^.(q) < 0.5) be- 
cause only a small fraction of the HI emission oc- 
curs on this scale (see Fig. [TT|) . Furthermore, the 
coefficients of the 70 /im-Ha correlation are higher 
than those of the 70 /im-neutral gas correlation on 
scales a < 6.3 kpc. 

6. Classical correlations between dust and 
gas 

The wavelet cross-correlations for different scales 
in Fig. [T^ show on which scales the distributions of 
the various types of emission are significantly cor- 
related. However, because the scale maps are nor- 
malized and information about absolute intensities 
is lost, they cannot be used to find quantitative re- 
lations between components of the ISM. Hence, to 
obtain numerical equations relating two distribu- 
tions, we need classical correlations. Classical cross- 
correlations contain all scales that exist in a distri- 
bution. For example, the high-intensity points of 
the Ha-70/im correlation in Fig. 1 141 represent high- 
emission peaks on small scales in the spiral arms 
(compare Fig. llOb ). whereas low- intensity points 
represent weak emission around and between the 
arms on larger scales (compare Fig. [TOjl) . The cor- 
relation coefficient of 79% is a mean of all scales, 
consistent with Fig. [T^l 

We made pixel-to-pixel correlations between the 
distributions of th^ and H2, HI, total gas as well 
as between de-reddened Ha and 24 /im, 70 /im, and 
160 /im. We restricted the comparisons to radii 
where all data sets are complete, R < 50' (or 
11.4 kpc), and to intensities above 2 x rms noise. 
To reduce the influence of the gradient in the gas- 
to-dust ratio (see Sect. 4.2), we calculated corre- 
lations for two radial ranges: 0' < i? < 30' and 
30' < R < 50'. We obtained sets of independent 
data points, i.e. a beam area overlap of < 5%, 
by choosing pixels spaced by more than 1.67 x the 
beamwidth. Since the correlated variables are not 
directly depending on each other , we fitted a powe r 
law to the bisector in each case (jlsobe et al.lllQQOl ). 

We also calculated the correlation coefficient, Tc, 
to show how well two components are correlated, 
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Fig. 10. Distribution of the de-reddened Ha emission (a) and the wavelet decomposition for scales 0.4, 
1.6, 4.0 kpc (6 to d). The central 2kpc was subtracted from the Ha map before the decomposition. The 
cross in the Ha map indicates the location of the center. 



and the student-t test to indicate the statistical sig- 
nificance of the fit. For a number of independent 
points of n> 100, the fit is significant at the 3cr 
level at > 3. Errors in intercept Uc and slope b of 
the bisector are standard deviations (Icr). 



We first discuss the correlations between the 
neutral gas and dust extinction, scaled from th„. 
Then we investigate the relationships between the 
emissions from dust and ionized gas. The results are 
given in Tables [5] and [Bl and examples of correlation 
plots are shown in Figs.[T51to[T5] 
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6.1. Correlation between neutral gas and dust 
extinction 

In search of a general relationship between neutral 
gas and dust extinction, a number of authors em- 



ployed scatter plots between gas column densities 
and extin ction, optical depth or FIR surface bright- 
ne ss (e.g. Savage et al . 1978; Boulangcr ct al. 1996; 
W alterbos fc Kennicutt. 1988; .Xu fc Hclou .1996 ; 
N eininger et al.lll998l:lNieten et al.ll2006[) . Thev ob- 
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Power-law fits to the various correlations given in 
Table 5. Thick lines: 30' < i? < 50'; thin lines: 
0' < i? < 30'. The shift between these intervals of 
the fits to the N(HI)-Ahq and N(gas)-AHa corre- 
lations is due to the radial increase in the atomic 
gas-to-dust ratio. 



tained nearly linear relationships between these 
quantities. As the studies on M 31 have various 
shortcomings (lower limits for extinction, H2 data 
not included and/or low angular resolution), we cal- 
culated classical correlations between the distribu- 
tion of dust extinction Ahq= 1.086 res and those 
of N(HI), N(2H2) and N(gas) at our resolution of 
45". As the correlations are restricted to gas column 
densities above 2x the rms noise, values of th„ < 
0.04 are not included (see upper panel of Fig. [T^ . 
The bisector fits given in Table 5 are plotted in the 
bottom panel of Fig. [T31 

The relationships between Ahq and N(2H2) for 
the two radial ranges are the same within errors, 



so the two areas can be combined. With a corre- 
lation coefficient of rc ~0.6, the correlation is not 
very good, indicating that only a small part of the 
extinction is caused by dust in molecular clouds. 
This is not surprising in view of the low molecular 
gas fraction in M 31 (see lower panels of Fig. [S]) 
and the small area filling factor of the molecular 
gas compared to that of the atomic gas. 

The correlations between Ana and N(HI) are in- 
deed better (rc ~ 0.7) than those between Ahq and 
N(2H2), but the relationships for the two radial in- 
tervals are not the same. Although both are nearly 
linear (power-law exponent b ~0.9), their power 
laws are shifted (see Fig. [T5)) in the sense that the 
values of Ahq in i? = 30' — 50' are about a factor 
of 2 lower than those inside R — 30'. This differ- 
ence is caused by the radial decrease of th„/N(HI) 
discussed in Sect. 4.2. The variation of this ratio 
within each of the radial intervals contributes to 
the spread in the scatter plots and reduces the cor- 
relation coefhcients. 

The correlations between total gas N(gas) and 
Ana are best (jc — 0.8), as Ahq represents dust 
mixed with both HI and H2 . They are close to linear 
{b ~0.8) and differ by nearly a factor 2 in Ana- 
The scatter plots for the two intervals are shown 
in the upper panel of Fig. 1131 In linear plots both 
power-law fits are going through zero, suggesting 
that the dust causing the extinction and neutral 
gas are mixed down to very low densities. 

Interestingly, extinction (or dust opacity) is pro- 
portional to the square root of N(2Il2), while it is 
about linearly related to the atomic gas density. 
This is due to the quadratic de pendence of N(2Ho ) 
on N(III) in M 31 observed by (jNieten et alll200gi . 
This dependence is expected if in cool, dense, and 
dusty HI clouds the formation and destruction rates 
of H2 are balanced (|Reach fc Boulangerlll998f) . 

6.2. Correlation between ionized gas and dust 

Because the emission from ionized gas is a good 
tracer of the present-day star formation rate and 
massive stars both heat the dust and ionize the gas, 
a correlation between the emissions from warm dust 
and ionized gas is expected. Relationships between 
the emission at 24 ^m and Paa or Ha emission from 
HII regions in nearby galaxies as well as relation- 
ships between glob al luminosities of gala xies have 
been reported (see iKennicutt et al.ll2009l and ref- 
erences therein). 

For M 31, the correlation between the emis- 
sio n from dust an d ioniz ed gas was first tested 
by iHoernes et al.l (|1998D . who found a good, 
nearly linear correlation between warm dust emis- 
sion and free-free radio emission for the radial 
range 30' < i? < 90', using HIRAS data and 
multi-wavelength radio data. Here we correlate 
the extinction-corrected Ha emission presented in 
Fig. [TUk with dust emission in the MIPS maps. 
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Table 5. Power-law relations and correlation coefficients Tc between dust extinction and gas components. 
Ordinary least-squares fits of bisector \og(Y)—ac + 61og(X) through n pairs of (logX, logY), where n is 
the number of independent points (Isobe et al. 1990); t is the student-t test. 
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The wavelet correlations in Fig. [12] (bottom- 
right) show that in M 31 Ha emission is best cor- 
related with dust emission at 70 /im. This suggests 
that of the MIPS bands, the 70 fim emission could 
best be used as the tracer of present-day star for- 
mation, making a numerical relation between the 
emissions at 70 fim and Ha of interest. Table |6] gives 
the bisector fits for the two radial ranges, which are 
very similar. Therefore, we present this correlation 
in Fig.[T4]for the entire radial range of 0' < i? < 50'. 
The power-law fit for this radial interval is 

log(/7o) = (-0.98 ± 0.02) + (1.09 ± 0.02) log(/H„), 

where /70 is in MJy/sr and /hq in 
10~^erg s~^ cni^'^ sr~^. The correlation is 
quite good (rc ~ 0.8) and nearly linear. In a linear 
plot the power-law fit goes through zero suggesting 
that the correlation is also valid for the lowest 
intensities. The good correlation indicates that the 
heating sources that power dust emission at 70 fim 
and ionize the gas must indeed be largely the same. 

Naturally, Ha emission is less correlated with 
the emission from cold dust at 160 /im than with 
emission from warm dust seen at the shorter wave- 
lengths. This is especially so at i? < 30' where the 
radial proffies differ most (see Fig. Moreover, 
the relation between the emission from cold dust 
and Ha is non-linear (see the bisector slope b in 
Table El) . 

Table \6\ shows that the correlations with 24 /im 
are slightly worse than those with 70 /im. In con- 
trast, in M 33 the 24/zm-Ha co rrelation is better 
than t he 70/Km-Ha correlation ([Tabatabaei et al.l 
l2007a[ ). This may suggest that in early- type galax- 
ies like M 31 the contribution from evolved AGB 
stars to the 24 /im emission is larger than in late- 
type galaxies like M 33. A significant stellar contri- 
bution to the 24 fim emission from M 31 is also in- 
dicated by the enhancement of the 24 ^m-to-70 fim 
intensity ratio in inter-arm regions where the radi- 
ation field is weak (Fig. [7]). 

Across M 31, the 24 /im emission is linearly 
proportional to the extinction-corrected Ha emis- 



sion (5 — 0.98 ± 0.02). A linear relationship was 
also found between the luminosities at 24/Km and 
extinction-corrected Paa of HII regions in M 51 
([Calzetti et al.l[2005[ ) and between the luminosities 
at 24 fim and extin ction-corrected Ha of HII regions 
in M 81 ([Perez-G onzalez et al, 2006). Comparing 
the 24 ^m luminosities and corrected Ha luminosi- 
ties of HII regions in 6 nearby g alaxie s (including 
M 51 and M 81), iRelano et all ([2007^ obtained a 
somewhat steeper power law with index 1.21 ±0.01, 
in agreement with the index for global lum inosities 
of galaxies (see also ICalzetti et al.l [2"007t) . Thus, 
while the L24 - Lhq relationship is linear within a 
single galaxy, the relationships for HII regions in a 
sample of galaxies and fo r global luminos it ies ar e 
non-hnear. According to iKennicutt et al.l ([20091 ) . 
the steepening is due to variations between galaxies 
in the contribution from evolved, non-ionizing stars 
to the heating of the dust that emits at 24 /Ltm. 

7. Star formation rate and efficiency 

Over the last 40 years many authors have stud- 
ied the relationship between the rate of star for- 
mation and gas density in M 31 by comparing the 
number surface density of mass ive you ng stars or 
of H II regions with that of HI (TBerkhuiisen 19771: 
Te njes fc HaudI [1991 lUnwinI 119801 : iNakai fc Sofud 
J98J). They found power-law exponents near 2 as 
was also obtai ned for the solar neighborhood by 
'Schmidt ([1959f ) , who firstj)roposed this relationship 
with HI volume density. [Kcnnicutt (1998a) showed 
that a similar relationship is expected between SFR 
and gas column densities. The early studies suffered 
from the effects of dust absorptio n and could not 
consi der molecular gas (apart from iTenies fc HaudI 
[T991I) . As the necessary data are now available, we 
again address this issue. 

We compared the distribution of the Ha emis- 
sion corrected for dust attenuation (see Fig. [TUb ) 
with those of HI, H2 and total gas. The corrected 
Ha emission is a good measure for the present-day 
star formation rate (SFR) that we first estimate 
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Table 6. Power-law relations and correlation coefficients between the emission from dust and ionized 
gas. Ordinary least-squares fits of bisector log(Y)=ac -I- 61og(X) through n pairs of (logX, logY), where n 
is the number of independent points (Isobe et al. 1990); t is the student-t test. 
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Fig. 14. Scatter plot between the surface bright- 
nesses of ionized gas and dust emission at 70 /um 
for the radial range 0'-50'. Only independent data 
points (separated by 1.67x beamwidth) with values 
above 2x rms noise are included. The line shows 
the power-law fit given in Table 6, which has an 
exponent close to 1. In a linear frame, this fit goes 
through the zero point of the plot. 



for the total area observed using the relation of 
(|Kennicuttlll998bD : 



SFR(M0yr-i) = 



iH„(ergs 1) 



(2) 



1.26 X 1041 ' 
where Lh„ is the Hq luminosity. In an area of 

llO.o' X 38.5' {R < 17kpc), the luminosity of 
the de-reddened Hq emission is in = 4.75 x 



10''°ergs-i or Lh„ = 1.7 x 10^ Lq for the dis- 
tance to M 31 of 780 kpc (see Table 1), giving 
SFR = 0.38 M0yr~i. However, this value is rather 
uncertain for two reasons. First, the contribution 
from the inner disk (i? < 25', nearly 6 kpc) is 
overestimated because in this area the number of 
ionizing stars is low and the gas must be mainly 
heated by other sources (see Sect. 7.1). Second, our 
Ha map is limited to about 55' (12.5 kpc) along 
the major axis, so some of the emission between 
i? = 12.5 kpc and R = 17 kpc is missing. Subtracting 
the luminosity from the area R <6 kpc gives a lower 
limit to the SFR of O.27M0yr-i for the radial 
range 6< R < 17 kpc. Earlier estimates of the re- 
cent SFR for a. larger part of the di s k ind i cated 0.35 
- lM ry.vr-1 (IWalterbos fc Braunl 11994 IWilliams ' 
120031: iBarmbv et all I2006D . Recentlv. iKang eTa! 
(|200l derived a SFR of 0.43 Mgyr"! (for metal- 
licity 2.5 X Solar) from UV observations of young 
star forming regions (<10Myr) within 120' from 
the center {R < 27 kpc). Their Fig. 13 suggests that 
about 20% of this SFR is coming from R > 17 kpc 
and a negligible amount from R <Q kpc. So for the 
range 6 kpc < R < 17 kpc they find a SFR of about 
0.34 MQyr~i, which is consistent with our lower 
limit of O.27M0yr~i. 

A SFR of 0.3 Mq yr~i yields a mean face-on sur- 
face density of Ssfr — 0.4 Mq Gyr"^ pc~^ between 
i? = 6kpc and i?=17kpc. This is about 6 times 
lowe r than the va l ue of S s f r = 2.3 Mq Gyr-ipc-2 
that IVerlev et all ()2009[ ) obtained for the disk of 
M 33 (i? < 7kpc), also using de- reddened Ha data. 

We can also calculate the star formation effi- 
ciency between it!=6kpc and i?=17kpc in M31. 
The total molecular gas mass in the entire area 
oi R < 17 kpc in M 31 is M(H2) = 3.6 x 10^ Mq 
(jNieten et al.l 120061 ) and that in the area 6 kpc < 
R <17kpc is M(H2) =2.9 X IO^Mq. Hence the 
star formation efficiency SFE = SFR/M(H2) be- 
tween 6 kpc and 17kpc radius is SFE = 0.9 Gyr"^. 
It is equivalent to a molecular depletion time scale 
of 1.1 Gyr. Hence, the disk of M 31 is about 
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three times less efficient in forming young 
massive stars than the northern part of the 
disk of M 33 (jGardan et al.l 120071 ). 

7.1. Star formation rate in the '10 kpc ring' 

The radial distributions of Ha emission in the 
bottom panels of Fig. [6] show a steep decrease 
from the center to i? ~ 2 kpc, followed by a shal- 
lower decrease to a minimum near i? = 6kpc. In 
the radial profile of Ssfr, shown in the upper 
panel of Fig. [T7] for the total area, the inner arms 
at i? = 2.5kpc and i? = 5.5kpc are only visible 
as little wiggles superimposed onto a high back- 
ground. Clearly, the starforming regions in these 
arms hardly c ontribute to the ioniza tion of the gas 
at R < 6kpc. iDevereux et al.l (|1994[ ) noted that at 
these radii the Ha emission is filamentary and un- 
like that in starforming regions, and since not many 
young, massive, ionizing stars arc found interior 
to th e '10 kpc ring' (Bcrkhuiiscn fc Humphreys! 
[T989t iTenies fc HaudI 119911: iKang et al.l I2OO90 the 
gas must be ionized by other sources. Naturally, 
the same holds for the heating of the warm 
dust, the emission of which also strongly increases 
towards the center. In an extensive discussion, 
IDevereux et al.l (|1994|) concluded that a collision 
with another galaxy in the past may explain the 
ionization of the gas and the heating of the dust 
as well as several other peculiarities (e.g. the dou- 
ble nucleus) in t he inner disk of M 31 (see also 
iBlock et al.ll2006[) . We note that the UV emission 
may also be influenced by this event, because it 
shows a similarly steep increase towards the center 
as the Ha emission. Furthermore, all radial profiles 
that increase towards the center are anti- correlated 
with the radial profiles of HI and total gas (see 
Fig. ini bottom panels). This leads to apparent de- 
viations in the Kennicutt- Schmidt law in the inner 
disk if EsFR is calct ilated from the usual star for- 
matio n tracers (see iBoissier et al.l [20071 : lYin et all 
[20091) . 

If massive stars are not responsible for the ion- 
ization of the gas and the heating of the dust, we 
can neither use the Ha emission nor the infrared 
emission as tracers of present-day star formation 
at -R < 6 kpc in M 31 , as was also pointed out by 
IDevereux et al.l ([1994!). Therefore, we investigated 
the relationship between SFR and neutral gas only 
for the interval 30' < i? < 50' (i?= 6.8-11.4kpc) 
containing the '10 kpc ring'. 

The correlation plots in Fig. [T5] and the results 
in Table [7| show that Ssfr is not well correlated 
with the surface densities of either H2, HI or to- 
tal gas (re ~ 0.45-0.59). In spite of this, the fit- 
ted bisectors are statistically significant {t > 3). 
Interestingly, we find a linear relationship between 
EsFR and S2H2 (exponent b = 0.96±0.03), which 
closely agrees with the average relationship for 7 
nearby galaxies, much brighter than M 31 (see 



Fig. [T5J1), analyzed bv iBigiel et all ^OM)- While 
in these galaxies molecular hydrogen is the domi- 
nant gas phase, most of the neutral gas in M 31 
is atomic (compare Fig. ITSk.b). Hence, the surface 
density of SFR is linearly related to that of molec- 
ular gas, irrespective of the fraction of molecular 
gas or the absolute val ue of the total gas surface 
density in a galaxy. Big iel et al.l (|200^ arrived at 
the same conclusion after comparing the galaxies in 
their sample. 

The correlation between Ssfr and total gas sur- 
face density is slightly better than that between 
SsFR and molecular gas surface density. The bisec- 
tor fit in Table 7 yields the Kennicutt-Schmidt-law 

SsFR = (0.076 ± 0.005) SJ;i°s±°■°^ (3) 

where Sqas and Ssfr are in Mq pc~^ and 
M0Gyr~^pc~^, respectively. The exponent of 
1.3 ± 0.05 is well in the range of 1.1-2.7 derived 
by iBigiel et"al] (|2008f ) . As a galaxy of low surface 
brightness, the SFRs in M 31 are correspondingly 
low. Our Esfr-Sqas relationship nicely fits on the 
low-brightness extension of the compilation o f avai l- 
able galaxy data in Fig. 15 of IBigiel etaD (|2008D . 
formed by the outer parts of their 7 galaxies and the 
global values for 20 galaxies of low surface bright- 
ness. 

Very recently, iBraun et al.l (|2009() also studied 
the dependence of SFR on gas density in M 31 us- 
ing th e new Westerbork H I survey and the CO sur- 
vey of lNieten et al] (|2006l ). They estimated the SFR 
from the surface brightnesses at IRAC 8 /im, MIPS 
2 4 Mm and GALEX F UV following the procedure 
of iThilker etHI (|2005f ). Our Fig.[Tia is comparable 
to the radial range 8-16 kpc in their Fig. 20D that 
shows the same range in Esfr as we find. No t e that 
the molecular gas densities of iBraun et al.l (|2009D 
are a factor of 1.6 larger (-1-0.21 dex) and have a 
larger dynamic range than our values due to dif- 
ferences in scaling of the CO data, inclination, an- 
gular resolution and radial ran ge. Scaling our rela - 
tionship to the assumptions of IBraun et al.l (|2009l) 
gives log(EsFR) = -0.44 + 0.961og(E2H2), which 
is in good agreement with their Fig. 20D. 

The dependencies of SFR surface density on to- 
tal gas surface dens ity in Fig. [T5h and in Fig. 20E of 
IBraun et all (|2009t ) have the same pear-like shape 
characterized by a broadening towards lower Esfr 
and a rather sharp cut-off near Eqas = 10 M0 pc~^. 
The cut-off comes from the Esfr-Ehi relation (see 
Fig. 15b) and occurs at t he same value as i n the 
bright galaxies analyzed bv lBigiel et al.l(|2008^ ■ who 
interpreted the lack of hi gher surfa ce mass densities 
as a saturation effect. Br aun et all (2009) show that 
in M 31 this truncation indeed vanishes after cor- 
recting the HI data for opacity, which could lead to 
somewhat steeper slopes in Figs. 15b and c. 
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Fig. 15. Scatter plots between the surface density of the star formation rate and neutral gas surface 
densities for the radial interval 30' < R < 50' (6.8 kpc < R < 11.4 kpc). All surface densities are face-on 
values. Only independent data points (separated by 1.67x beamwidth) above 2x rms noise were used, 
a. SFR versus molecular gas; b. SFR versus atomic gas; c. SFR versus total gas. Full lines indicate 
the power-law fits given in Table 6. The relationship in a. is linear and nearl y the same as the average 
relationship for 7 bright galaxies (dashed line) derived bv iBigiel et al.l (|2008f ). Note that the cut off in 
Ehi and Egas is near 10 Mq pc~^. 



Table 7. Kennicutt- Schmidt law in M 31 for 30' < R < 50'. Ordinary-least squares fits of the bi- 
sector log(I]sFR)=ac + &log(E), where Ssfr is the face-on surface density of the star formation rate 
in Mq Gyr~^ pc~^ and E the face-on value of the gas surface density in Mqpc"^; n is the number of 
independent points; Tc is the correlation coefficient and t the student-t test. 
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7.2. Radial variations in the Kennicutt-Schmidt law 

In Fig. ITBk . we plot the mean values in 0.5 kpc- 
wide rings in the plane of M 31 of Esfr against 
those of Egas from i? = 6kpc to i?=16kpc. 
The points form a big loop with a horizontal branch 
for i? = 6 - 8.5 kpc and a maximum Esfr in the 
ring i?=10.5 - 11.0 kpc (see also Fig. [13 This 
beh avior was already note d by iBerkhuiisenI (|l977t ) 
and iTenies fc HaudI (|1991[ ) , who used the number 
density of HIT regions as tracer of SFR and HI 
gas, a nd was recently confirmed by iBoissier et al.l 
(I2007D from GALEX UV data and total gas 
Both lBerkhuiisenI (|1977D and lTenies fc HaudI (|l99lf ) 
showed that the differences between the slopes in- 
side and outside the maximum of the starforming 
ring is greatly reduced when the increase in the 
scale height of the gas with increasing radius is 
taken into account. We calculated the scale height, 
h, from th e scale height of the HI gas given by 
Eq. 13 of iBraunI (|1991[ ). scaled to D = 780 kpc, 
assumed half this value for that of the H2 gas, 



and a constant scale height for the ionizing stars. 
Fig. 116b shows Esfr as a function of gas volume 
density riGAS = N(HI)/2h + N(2H2)/h. The points 
have moved towards each other, but the horizontal 
branch remained and the behavior on the starform- 
ing ring has become more complicated. 

The variations in slope in Fig. [TB] are clear ev- 
idence for radial variations in the index of the 
Kennicutt-Schmidt law. Such variations are not 
specific to M 31 as th ey are also se e n in s ome of the 
galaxies analyzed bv iBigiel et al.l (|2008[ ). In order 
to quantify the variations, we determined the bisec- 
tors in scatter plots for three circular rings: R = 7 
- 9 kpc, R — 9- 11 kpc and _R=11 - 13 kpc, cover- 
ing the horizontal branch, the increasing part in- 
side the maximum and the decreasing part outside 
the maximum, respectively. The results are given in 
Table 8. The index for the star fromation law for 
surface densities is unity for the 7-9 kpc ring and 
about 1.6 for the other two rings. Hence, the slope 
of &= 1.30 ±0.05 obtained for the '10 kpc ring' 
(i? = 6.8- 11.4kpc) in Sect. 7.1 represents the mean 
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value of the first two rings considered here. The 
scatter plots between Ssfr and gas volume den- 
sity yield bisector slopes that are about 0.2 smaller 
than those for surface density. The correlation co- 
efficients are all close to Vc = 0.59 ± 0.03 for the 
'lOkpc ring' (see Table 7), indicating that even 
in 2 kpc-wide rings the intrinsic scatter is consid- 
erable. This implies that on scales of a few hun- 
dred parsec significant variations in the index of 
the Kennicutt-Schmidt law and in the star forma- 
tion efficiency occur. 

7.3. Radial variations of SFR and SFE 

In Fig. [T7] (upper panel) we present the radial pro- 
file of the SFR surface density between 6kpc and 
17kpc, averaged in 0.5 kpc-wide circular rings in 
the plane of M 31. The face-on value s vary between 
about 0.1 and 1 Mf^ Gyr~^ pc~^. iBoissier et al.l 
(|2007[) and lBraun etai] (f2009[ ) obtained similar val- 
ues for Ssfr in this radial range from GALEX UV 
data. They are about 10 times smaller than the sur- 
face densities of SFR between R= 1.5 kpc and 
R = 7kpc in the north ern part of M 33 observed 
by (|Gardan et al.ll2007l) . 

In the lower panel of Fig. [17] we show the ra- 
dial profiles of the surface density of the molec- 
ular gas and of the star formation efficiencies 
SFE = Esfr/E2H2 and Esfr/^gas- Although the 
maximum Esfr, occurs on a relative maximum in 
the molecular gas density (in ring 10.5-ll.Okpc), 
Ssfr is only about 70% of its maximum value 
where the molecular gas density is highest (in ring 
9.0-9.5 kpc). Consequently, SFE varies significantly 
with radius. Between R — 6 kpc and R — 15 kpc 
SFE fluctuates around a value of 0.9 Gyr^^, with 
a minimum of 0.46 ± 0.01 Gyr~^ near i?=9kpc. 
Thus SFE is smallest where E2H2 is highest! Up 
to i?=12kpc the efficiency Ssfr/Sqas shows the 
same trend as SFE. The increase in SFE between 
12 kpc and 15 kpc radius of a factor 1.5 results 
from the difference in radial scale lengths of Esfr 
(or Ha emission) and the molecular gas density 
(see Fig. [6l and Ta ble 3). Interestingly, in M 33 
([Gardan et al.|[2007l ) found a radial increase in SFE 
of a factor 2 between 2 kpc and 6 kpc radius with 
similar fluctuations around the mean as we observe 
in M 31, but the mean value in M 31 is abo ut three 
times lower than in M 33. Furthermore, Ler ov et al.l 
( 2008h found signiflcant variations in the efficiency 
Esfr/Sqas on a linear scale of 800 pc in the sample 
of 12 spiral galaxies analyzed by them. 

That large, small-scale variations in SFE exist 
in galaxies is also clear from the large spread in the 
scatter plots of I]sf r.~S2Hi visib l e in F ig. [T5k and 
in several figures of lBigiel et all (|2008h . The same 
value of Ssfr can occur in a range of S2H2 spanning 
more than a factor of 10. 

We may conclude that neither the present-day 
star formation rate Ssfr nor the star formation effi- 
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Fig. 17. Radial variation of the face-on surface den- 
sity of the star formation rate Ssfr and star for- 
mation efficiency SFE in M 31, averaged in 0.5 kpc- 
wide rings in the plane of the galaxy. Beyond 
R = 12 kpc the data are not complete, because the 
observed area is limited along the major axis (see 
Fig. 10a). Top: Radial profile of Ssfr- Bottom: 
full line - SFE = SSFR/S2H2; long dashed line - 
Ssfr/Sqas; dots - radial profile of the molecular 
gas surface density S2H2 seen face-on. Note that 
Ssfr and S2H2 do not peak at the same radius. 
Statistical errors in Ssfr and S2H2 are smaller than 
the symbols and those in SFE and Ssfr/Sqas are 
smaller than the thickness of the lines. Only be- 
yond i? = 12 kpc the error in SFE slowly increases 
to 0.3 Gyr"^ at i?=15kpc. 



ciency SFE is well correlated with the molecular gas 
surface density. Hence, other factors than molec- 
ular gas density must pla y an important rol e in 
the star formation process. iBigiel et all ()2008f ) ar- 
gue that local environmental circumstances largely 
determine the SFE in spiral gal axies. These factor s 
are extensively discussed by e.g. lLerov et al.l (|2008[ ). 
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Fig. 16. Mean face-on values of Ssfr, averaged in 0.5 kpc-wide circular rings in the plane of M 31, plotted 
against the corresponding mean values of (a) gas surface density SqaSj and (b) gas volume density ncAS- 
The upper left point is for the ring i? = 6.0-6.5 kpc, the maximum in Ssfr is in ring 10.5-11.0 kpc and 
the minimum in ring i? = 14.5-15.0kpc. Typical errors are 0.01 Mq Gyr~^ pc~^ in Ssfr, O.O2M0pc~^ 
in SqaS: and 3 x 10~^ Mq pc~^ in riQAS- Points for R > 12kpc suffer from missing data points near the 
major axis, the number of which increases with radius. 



Table 8. Kennicutt-Schmidt law in three radial intervals in M 31. Ordinary least-squares fits of the 
bisector log(I]sFR)=ac + ^log(X), where Ssfr is the face-on surface density of the star formation rate in 
M0 Gyr^^ pc~^, X = Eqas is the face-on gas surface density in M0 pc~^ and X = tiqas is the gas volume 
density in at cm~'^. n is the number of independent points; Tc is the correlation coefficient and t the 
student-t test. 



i?(kpc) 


X 




b 


n 


rc 


t 


7-9 


Sgas 


-0.90±0.04 


1.03±0.06 


216 


0.54±0.06 


9 


9-11 




-1.43±0.06 


1.67±0.08 


356 


0.63±0.04 


15 


11-13 




-1.46±0.06 


1.55±0.08 


297 


0.62±0.05 


13 


7-9 


wgas 


0.18±0.04 


0.88± 0.06 


218 


0.51±0.06 


9 


9-11 




0.45±0.04 


1.50±0.07 


356 


0.62±0.04 


15 


11-13 




0.35±0.04 


1.35±0.07 


297 


0.62±0.05 


14 



8. Summary 

In this paper, we studied the emission from dust, 
neutral gas and ionized gas in the disk of M 31, and 
the relationships between these components on var- 
ious linear scales. We compared the Spitz er MIPS 
maps at 24 /im, 70 /xm and 160 /zm ( ,Gordon et al.l 
|2006[) to the distributions of ato mic gas seen in the 
HI line (jBrinks fc Shand Il984l), molecular gas as 
traced by the CO (1-0) line (iNieten et al.l l2006D 
and ionized gas observed in Ha ( Devereux et al.l 
1199411 . All data were smoothed to an angular res- 
olution of 45" corresponding to 170 pc x 660 pc in 
the plane of the galaxy. 

For each of the dust and gas maps, we calcu- 
lated the mean intensity distribution function 
of radius (Fig. IB]), separately for the northern and 
the southern half of M 31. Using wavelet analysis, 



we decomposed the dust and gas distributions in 
spatial scales and calculated cross-correlations as a 
function of scale. We also used classical correlations 
to derive quantitative relations between the various 
dust and gas components. 

Using the MIPS 70 /im and 160 /xm maps, we 
derived the distributions of the dust temperature 
and optical depth. The dust optical depth at 
the Ha wavelength was used to a) investigate 
the dust-to-gas ratio, b) derive scaling relations 
between extinction and neutral gas emission, and 
c) de-redden the Ha emission in order to estimate 
the recent star formation rate. We also presented 
the Kennicutt-Schmidt law indices obtained for 
the bright emission ring near R— 10 kpc in M 31. 
We summarize the main results and conclusions as 
follows. 
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1. Dust temperature and opacity: 

• The dust temperature steeply drops from about 
30 K in the center to about 19 K near R = 4.5 kpc, 
and stays between about 17 K and 20 K beyond 
this radius (Fig. 2). The mean dust tempera- 
ture in the area studied is about 18.5 K. This 
is 3 K less than the tempe rature obtained by 
IWalterbos fc SchweriiiS ^98?) between the IRAS 
maps at 60 /xm and 100 ^m that both trace warmer 
dust than the MIPS maps at 70 /xm and 160 /xm 
used here. 

• The dust optical depth at Ha along the line of 
sight varies in a range between about 0.2 near the 
center and about 1 in the '10 kpc ring' (Fig. 4) 
with a mean value of 0.7±0.4 (the error is standard 
deviation) and a most probable value of ~0.5, 
indicating that M 31 is mostly optically thin to 
the Ha emission. The total flux density of the 
Ha emission increases by 30% after correction for 
extinction. 

2. Radial distributions: 

• The radial scale lengths between the maximum 
in the '10 kpc ring' and i? = 15 kpc of the warm 
dust are smaller than that of the cold dust, as is 
expected if the warm dust is mainly heated by UV 
photons from star forming regions and cold dust 
by the ISRF. With the largest scale length, atomic 
gas has the largest radial extent of the dust and 
gas components considered here. 

• The radial gradient of the total gas-to-dust ratio 
is consistent with that of the oxygen abundance 
in M 31. The gas-to- dust ratios observed in the 
solar neighborhood (|Bohlin et al.l I1978D occur 
near R = 8.5 kpc in the disk of M 31 where 
N(gas)/rHa= 2.6 x 10^^ atcm"^. 

3. Properties as a function of scale: 

• Spatial scales larger than about 8 kpc contain 
most of the emitted power from the cold dust and 
the atomic gas, whereas the emissions from warm 
dust, molecular gas and ionized gas are dominated 
by scales near 1 kpc, typical for complexes of star 
forming regions and molecular clouds in spiral 
arms (Fig. 11). 

• Dust emission is correlated (r^, > 0.6) with both 
neutral and ionized gas on scales > 1 kpc. 

• On scales < 1 kpc, ionized gas is best correlated 
with warm dust and neutral gas (both HI and H2) 
with cold dust. On the smallest scale of 0.4 kpc, an 
Hl-warm dust correlation hardly exists (r^, c± 0.4) 
because not much HI occurs on the scale of star 
forming regions (see Fig. 11). 

4. Relationships between gas and dust: 

• Ha emission is slightly better correlated with 
the emission at 70 /um than at 24 /im (Fig. 13, 
Table 6), especially on scales <2 kpc (Fig. 12). 
As in M 33 the 24/im~Ha correlation is best, this 



suggests that in early- type galaxies like M 31 the 
contribution from evolved AGB stars to the 24 /j,m 
emission is larger than in late-type galaxies like 
M 33. 

• Dust extinction Ana is not well correlated with 
N(2H2) indicating that dust mixed with molecular 
clouds does not contribute much to the total 
extinction. Although the correlation with N(HI) is 
better, Aho is best correlated with N(HI-|-2H2). 

• Dust opacity is proportional to the square root 
of N(2H2) but about linearly rela t ed to N(HI), 
as was also found by iNieten et all (|2006f) at 90" 
resolution. This is an indirect indication of a 
balance between the formation and destruction 
rates of H2 in cool, dusty HI clouds. 

• In the central 2 kpc both the dust opacity and 
the HI column density are very low and the dust 
temperature is high. This combination may explain 
the lack of H2 in this region. 

5. SFR and SFE: 

• The SFR in M 31 is low. The total SFR in the 
observed field between R — 6 kpc and i? = 17 kpc 
is O.27M0yr~^ and the star formation efficiency 
is 0.9Gyr~^, yielding a molecular depletion time 
scale of 1.1 Gyr. This is about three times longer 
than observed in the northern part of M 33 
(Gardan et al. 2007). The radial distribution of 
5^SFR in 0.5 kpc- wide rings in the plane of the 
galaxy (Fig. 17) varies between about 0.1 and 
IM0 Gyr~^ pc~^, values that are about 10 
times smaller t han i n the northern part of M 33 
(jGardan et al.1 I2007D . Between R — 6 kpc and 
i? = 15 kpc, SFE varies between about 0.5Gyr~^ 
and 1.5Gyr~^, whereas the efficiency with respect 
to the total gas surface density slowly decreases 
from about 0.18 Gyr'^ to about 0.03 Gyr"^ 

• SFR is not well correlated with neutral gas and 
worst of all with molecular gas in the radial range 
30' — 50' containing the '10 kpc ring' (Fig. 15, 
Table 7). In spite of this, the power-law fits are sta- 
tistically significant. We find a linear relationship 
between the surface densities of SFR and molecular 
gas (power-law exponent 0.96 ± 0.03), and a power 
law with index 1.30 ± 0.05 between the surface 
densities of SFR and total gas. These results 
agree with the average relationship f or 7 nearby 
galaxies much brighter than M 31 (jBigiel et al.l 
2008). While in these galaxies molecular hydrogen 
is the dominant gas phase, most of the neutral 
gas in M 31 is atomic. Thus, the surface density 
of SFR depends linearly on that of molecular gas 
irrespective of the fraction of molecular gas or the 
absolute value of the total gas surface density in a 
galaxy. 

Some important implications of this study are: 

- Precaution is required in using the total IR lu- 
minosity (TIR) as an indicator of recent SFR or 
to derive dust opacity for an early-type galaxy 
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like M 31, because the cold dust is mainly heated 
by the ISRF and the warm dust emission at 
24 ^m is partly due to evolved stars (especially 
in the bulge of the galaxy). 
- Neither the present-day SFR nor SFE is well 
correlated with the surface density of molecular 
gas or total gas. Therefore, other factors than 
gas density must play an important role in the 
process of star formation in M 31. 
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